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Electron reflection at an interface is a fundamental quantum transport phenomenon. The most famous 
electron reflection is the electron^hole Andreev reflection (AR) at a metal/superconductor interface. While 
AR can be either specular or retro- type, electron-^ electron reflection is limited to only the specular type. 
Here we show that electrons can undergo retro-reflection in bilayer graphene (BLG). The underlying 
mechanism for this previously unknown process is the anisotropic constant energy band contour of BLG. 
The electron group velocity is fully reversed upon reflection, causing electrons to be retro-reflected. Utilizing 
a BLG/superconductor junction (BLG/S) as a model structure, we show that the unique low energy 
quasiparticle nature of BLG results in two striking features: (1) AR is completely absent, making BLG/S 
100% electron reflective; (2) electrons are valley- selectively focused upon retro-reflection. Our results 
suggest that BLG/S is a valley-selective Veselago electron focusing mirror which can be useful in valleytronic 
applications. 

Specular reflection of a particle by smooth solid wall is a classic textbook illustration on the principle of 
conservation of momentum^ Because the particle motion is intrinsically 'locked' to the direction of the 
momentum, the conservation of momentum parallel to the interface, ky, automatically implies the group 
velocity component parallel to the interface, Vy, to remain unchanged upon reflection. This results in specular 
reflection (Fig. la). In the quantum mechanical case, retro reflection can occur at a superconducting interface via 
Andreev reflection (AR)^. AR is a two particle process in which an incident electron in the normal metal couples 
with another electron below the Fermi level to form a Cooper pair, crossing the interface into the superconductor^ 
and leaving behind a hole. Since the hole is a time reversed version of the incident electron, its Vy is anti-parallel 
with ky. The conservation of ky requires a sign reversal of Vy upon reflection. This results in electron^hole retro 
reflection. In graphene^ and a 2D semiconductor with Rashba spin-orbit coupling^ the pseudo-spin (graphene) 
and real-spin (Rashba semiconductor) nature of the two systems creates an additional 'spin'-split branch: up- 
down Dirac cones in graphene and left-right sub -bands in the Rashba semiconductor. When forming an interface 
with a superconductor, the split branches form two 'flavors' of holes: one with ky parallel with Vy and one with ky 
anti-parallel with Vy. This allows both retro and specular AR to occur in these systems. 

Retro reflection of electrons at an interface remains elusive. For retro electron reflection (RER) to occur, 
reversal of Vy is required but this, however, violates the momentum conservation principle. To achieve RER, Vy 
needs to be 'unlocked' from ky such that the direction of Vy can be freely reversed upon reflection while conserving 
ky (Fig. lb). We define a constant energy hand contour as the outline of the energy spectrum across a constant 
energy slice in /c-space. The concept of constant energy band contours will be useful in the qualitative under- 
standing of the electron reflection trajectory, since it is directly related to the electron group velocity. Fig Ic shows 
the circular band contour of a parabolic (or linear) energy spectrum. In this case, specular electron reflection 
(SER) is the only permissible electrons electron reflection process since Vy does not change sign upon reflection. 
We articulate that (i) anisotropic and (ii) opposite band contour for incident and reflected electron states are the 
fundamental criteria for an electron to be retro -reflected. Criterion (i) allows the electron motion to be unlocked 
from the direction of momentum since in an anisotropic band contour, the momentum no longer aligns with the 
group velocity while the direction of the group velocity is solely determined by the contour of the band via 
v= {vx, Vy) =h~^Vk£k' Criterion (ii) allows the direction of Vy to be reversed upon reflection. In Fig. Id, a 
hypothetical energy spectrum with boomerang- shaped constant energy band contour in A;-space is shown. 
The boomerang-shaped constant energy band contour is highly anisotropic in /c-space. The constant energy 
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Figure 1 | (a) Since electron group velocity component Vy is locked to the momentum component ky, conservation of ky allows only specular electron 
reflection (SER) to occur; (b) retro electron reflection (RER) requires y^to be 'unlocked' from ic^such that reversal of y^can still occur while conserving 
icp (c) constant energy slice of a parabolic (or linear) energy spectrum. The circular band contour allows only SER to occur; (d) constant energy slice of a 
hypothetical energy spectrum with boomerang-shaped band contour. The constant energy band contour of the incident states in /c-space is convex while 
that of the reflected states is concave. The opposite band contours between incident and reflected states causes the sign reversal of upon reflection. This 
results in RER. 



band contour of the incident states is convex while that of the 
reflected states is concave. Along a constant line of ky^ the reflected 
electron motion is completely reversed in both x- and y- direction. Vy 
is reversed whfle conserving ky. The reflection process is hence RER. 

The boomerang- shaped constant energy band contour occurs in a 
Bernal-stacked bilayer graphene (BLG). The low energy spectrum of 
BLG exhibits the trigonal warping effect (TW), i.e. the energy spec- 
trum splits into four discrete low energy pockets which join to form a 
single-band at higher energy^. The TW has resulted in highly aniso- 
tropic constant energy band contour which fulfifls both of the RER 
criteria. In the tight-binding framework, the low energy effective 
Hamiltonian is given as 



0 kl 
k\ 0 



0 



0 



where = 1(- 1) for K (K') vafley, Sq = {v^lvpfyi. k± = k^± iky. The 
wave vectors are normalized by ko and are dimensionless, i.e. k^ = kj 
koy ky = ky/ko where ko = v^yi/vp. The hopping parameters are yi = 
0.3 eV, Vs/vp = 0.1 and Vp = 10^ m/s^. The characteristic wave vector 
ko defines the phase- space separation between central and satellite 
Dirac points. The low energy pockets join to form a single trigonally 
warped band structure at Sk = So/4. The quadratic term in Eq. (1) 
represents the indirect electron hopping between B| and A| sublat- 
tices (t and | denote upper and lower layer respectively) mediated by 
A|-B| dimer pair whfle the linear term represents the weak direct 
electron hopping between the B| and A| side (Fig. 2a). The inclusion 
of the linear term anisotropically deforms the low energy dispersion: 
^k,^ = So \/k'^ + k^ — 2^k^ cos 3(/> where (j) is the momentum angle and 
= ± is the vafley- index. For simplicity, we denote Sk = 8k,+ for the 
K- vafley. The BLG energy spectrum exhibits three distinct categories 
of band contour (Fig. 2b): (A) the band structure splits into four 
distinct low energy pockets' at 8^ < eo/4; (B) the pockets join to form 
a single band with 'three-leaf -like contour at 8o/4 < Sjt < 10.9£o; and 
(C) the band contour becomes parabolic-like in the high energy 
regime of 8^ > 10.9£o- Regime (B) is the most interesting since the 
reflection process is dominantly RER. The 'three-leaf '-shaped band 



contour is highly anisotropic. The incident electron states residing on 
the right-pointing 'leaf and the right-edges of the upper and lower 
leaves has a convex constant energy band contour while the reflected 
electron residing on the left-hand edge of the upper and lower leaves 
has a concave constant energy band contour. As shown in Fig. 2c, an 
incident electron with direction of motion pointing upwardly-right 
(red arrow) is reflected as an electron pointing downwardly-left 
(green arrow) along a constant line of ky. The reversal of Vy whfle 
conserving ky results in RER. SER occurs only at large \ky\ in which 
criteria (ii) no longer holds (incident and reflected electrons has the 
same constant energy band contour and hence Vy does not change 
sign upon reflection). The very low energy regime (A) and high 
energy regime (C) are less ideal for RER. In Regime (A), RER occurs 
only minimally in the upper and lower pockets. In regime (C), the 
quadratic term in Eq. (1) dominates and the constant energy band 
contour loses the required boomerang-shaped constant energy band 
contour. 

In this work, we utflize a hybrid junction made up of a bilayer 
graphene and a superconductor (BLG/S) to theoretically dem- 
onstrate the unusual RER phenomena. The RER in BLG/S represents 
one last missing piece in the quantum transport phenomena of a 
superconducting interface where specular reflection of electrons, 
and retro and specular reflection of holes have afl been demon- 
strated^'^'^"^\ Not only RER occurs in such a junction, there are 
two additional surprises: (1) the 2n Berry phase of the BLG low 
energy electron forbids its transmission across the junction, resulting 
in the complete absence of AR; and (2) electron reflection in BLG/S 
is vafley- dependent. The BLG/S acts as a valley-selective Veselago 
electron-mirror. 



Results 

Determination of electron reflection type using group velocities. 

In the usual cases (parabohc or linear band structure), quasi- 
momentum angle (j) = tan"^ {kylk^ is commonly used to described 
particle trajectory. Each uniquely defines a particle traveling in the 
direction of In BLG, ^ is no longer locked to the electron motion 
due to the constant energy band contour anisotropy. Instead, the 
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Figure 2 | (a) Bernal-stacked bilayer graphene lattice structure; (b) the energy spectrum contour plot in phase space, showing three distinct energy 
regime: Regime (A): Sj^ < sJA, Regime (B): sJA^ < and Regime (C): > 10.9ifo. Retro electron reflection occurs optimally in Regime (B) due 

to its 'boomerang-like' anisotropic constant energy band contour in space. In high energy Regime (C), retro reflection is no longer possible as the band 
contour becomes parabolic-like; (c) Band contour of an constant energy slice in Regime (B). The green and blue arrows denotes incident and reflected 
electron direction of motion respectively. The opposite band contour between incident (convex) and reflected states (concave) causes Vy to reverse its 
direction upon reflection, leading to RER. 



electron trajectory is more accurately described by group velocity 
angle 0 = tan"^ (Vy/v^^)- Unlike (j), the uniqueness of 0 is lost. Each 
0 can simultaneously represent several electron states and hence 
becomes ambiguous. We therefore characterize the reflection 
problem by momentum component ky instead of the incident 
quasi-momentum angle 0 or incident group velocity angle (j). The 

group velocity is given as Vv = /z~ ^ dsj^'^ j dk^ where v = x,;/ denotes x- 
and )/- directional group velocity respectively. We consider the case 
where the reflection interface is placed along the y- direction. In 
Fig. 3a-c, Vx of incident and reflected electrons across a constant 
energy slice are shown, has the opposite sign between incident 
electron and reflected electron and is highly asymmetric due to the 
broken left-right symmetry of the BEG band contour. 

The group velocity component Vy is the most important quantity 
in determining the electron reflection type (Fig. 3d-f). Sign reversal 
of Vy at a fixed ky indicates RER (enclosed by dotted lines in the 
Fig. 3d-i). The incident electron Vy reaches a sharp maximum when 



the band contour is flattened as ky increases. At ky > 1.3/co, Vy no 
longer changes sign upon reflection and the reflection becomes 
specular. At small 8^, almost all permissible reflection states are 
within the RER window and hence RER is dominant. The specular 
reflection states (outside the RER window) expands as 8^ increases, 
making RER relatively less profound. At 8k = Ssq (Fig- 3f), more than 
half of the permissible ky corresponds to specular reflection. 

The non-uniqueness of 6 is obvious in Fig. 3g-i, i.e. multiple sets of 
Vx and Vy produce the same 9. In contrast to Vy, RER requires the sign 
of 0 to remain unchanged upon reflection since in an retro reflection 
the electron is reflected to the same side of the normal of incidence. 
Due to the asymmetric constant energy band contour in /c- space, the 
retro-reflected electrons do not trace back the incident path. Instead, 
they are mostly retro-reflected with smaller angle, maximally about 
20°. The reflection angle approaches zero with increasing 8^ as the 
constant energy band contour becomes 'smoother' and loses the 
anisotropic band contour required to effectively produce large-angle 
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Figure 3 | Group velocities of incident and reflected electrons, (a)-(c) Group velocity (d)-(f) group velocity v^; and (g)-(i) incident and reflection 
angles of an incident K valley electron at a given ky. The energies are ejt = 0.5eo> = 2eo and = 5eo respectively for column 1, 2 and 3. In (a)-(c), is 
positive for incident electrons and negative for reflected electrons. The RER regime is enclosed between the dashed lines in (d)-(i). In (d)-(f), sign reversal 
of corresponds to RER. In (g)-(i), the reflection angle of RER does not change sign since electron is reflected to the same side of the normal. At small 
almost all of the permissible reflections are RER. As E]^ increases, more permissible states becomes specular reflection states (which lies outside the RER 
windows) and the RER angle approaches 0°. 
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Figure 4 | RER in BLG/S junction. Green, blue and gray arrows indicate 
incident electron, retro-reflected electron and transmitted quasiparticle 
respectively. Transmission across the junction is forbidden due to the 2n 
Berry phase nature of bilayer graphene electron. The junction is hence 
100% electron-reflective. 

RER. At 8k > 10.9eo, the reflection 0 drops below zero, i.e. reflected 
electrons are placed on the opposite side of the incident normal, and 
therefore all electron reflection becomes specular-type. 

Absence of Andreev reflection at BLG/S. A highly electron- 
reflective interface is preferred for the observation of RER. Such an 
interface conveniently occurs in a ballistic BLG/S junction (Fig. 4) 
where electron transmission across the junction is completely sup- 
pressed due to the spinor structure of the BLG electron wave- 
function. The superconductor is assumed to be s-wave type with 
isotropic pairing potential A* = A and is induced intrinsically in 
BLG via the proximity effect^^ (thus avoiding lattice mismatch 
problems). The quantum transport across the junction can be 
described by Bogoliubov-de Genne (BdG) equations in the form of 
H^{r) = £^(r), where H is the Hamiltonian describing the 
excitation in BLG/S junction. The reflection coefficients can then 
be obtained by solving the BdG equation and matching the wave- 
functions at the boundary^^ (see Supplementary Information). The 
reflection probabilities were found to he Rh = 0 and Re = I for AR 
and electron-^ electron reflection respectively. Electron-^ electron 
reflection is therefore the only permissible process in the low ener- 
gy regime of BLG/S and such a junction acts as a 100% reflective 
electron-mirror. 

BLG/S junction as a valley-selective Veselago electron-mirror. The 

nature of RER in BLG/S is related to Veselago optics A Veselago 
electron-lens based on single layer graphene (SLG) and BLG p-n 
junctions has previously been reported^^'^^. Such junctions are 
Veselago electron-lenses because of the negative electron refractive 
index. In SLG, the negative electron refractive index originates from 
the existence of a lower Dirac sub -band which has opposite qua- 
siparticle dynamics in comparison with the upper Dirac sub-band^^. 
In BLG, the negative electron refractive index originates from its 
anisotropic constant energy band contour^^. In contrast to a Veselago 
electron- /ens, BLG/S is equivalent to a Veselago electron- m/rror. 
In optics, ray refraction is governed by Snell's Law 



Hi sin Oi = — fit sin 0^ 



(2) 



where rii and Ut are the refractive index at the incident side and at 
the transmitted side respectively. A negative sign is added to the 
right hand side of Eq. (2) since we define the reflection angle Oy to 
have the same sign as the incident angle Oi if the reflected ray is on 
the same side of the normal of incidence as the incident ray. For 
simplicity, we consider a Veselago lens with Uf = —fii. Since Oi = 
Or the transmitted ray is placed at the wrong' side of the normal 
of incidence and this results in the transmitted ray, originating 
from a point source, being re-focused after propagating a certain 
distance in the transmission side of the interface (Fig. 5a). The 
Veselago lens is therefore a flat lens with ray-focusing behavior. 
For the reflection problem, a reflection 'Snell's Law' can also be 




Figure 5 | Schematic drawings of (a) Veselago lens with /if = — Uy; and (b) 
Veselago mirror n^. ~ — /jy. In the Veselago lens, a ray emitted from a point 
source (denoted by yellow triangle) is focused at the transmitted side of the 
interface. For the Veselago mirror, ~ — tii is chosen for better visual 
clarity. The retro -reflected ray is focused at the incident side of the interface. 

written as Ui sin Of = —rir sin Or where rir is the refractive index on 
the reflected side. The reflection 'Snell's Law' is however redund- 
ant in the usual case, since reflection occurs in the same medium 
where Ur = tii and hence Or = —Of. To construct a Veselago 
mirror, we consider a hypothetical anisotropic medium with Ur 
= —Hi, i.e. the refractive index of the medium undergoes a sign 
change when a ray travels in the opposite direction with respect to 
the incident ray. This leads to the retro -reflection condition of Or 
= Of. The reflection interface is Veselago-like since the reflected 
ray is placed at the 'wrong' side of the incident normal. The retro - 
reflected rays, emitted from a point source, trace back their ori- 
ginal incident path and are re-focused at the incident side of the 
interface (Fig. 5b). The interface can therefore be regarded as a flat 
Veselago mirror with ray-focusing behavior. It becomes obvious 
that the RER effect in BLG/S is equivalent to a Veselago electron - 
mirror. In BLG/S, the absence of AR makes the interface a perfect 
electron -mirror and the RER process is Veselago-like. BLG/S is 
therefore a flat Veselago electron -mirror with electron focusing 
ability. 

In BLG, the low energy electrons reside in two inequivalent K and 
K' valleys. Intervalley scattering is strongly suppressed^^. Devices util- 
izing this robust valley degree of freedom in graphene, valleytronics, 
have been proposed^^. The first step towards valleytronics is a valley 
polarizer. The RER phenomena of BLG/S junction offers another 
simple way to separate electrons from different valleys without the 
need of cutting^^, straining^^, terahertz laser-irradiating^^ or selective- 
damaging^^ graphene. The K and K' electron excitation spectrum of 
BLG/S has an opposite constant energy band contour (Fig. 6a and 
Fig. 6b). The reflection states in the iC valley have a much 'smoother' 
constant energy band contour (Fig. 6a) than the reflection states in the 
K' valley (Fig. 6b). This immediately suggest that the electron 
reflection in BLG/S is valley-selective. In Fig. 6c, it can be seen that 
the iC electrons are predominantly retro -reflected with small angle 0 < 
20° while the K' electrons are predominantly retro-reflected with large 
angle 50° < ^ < 70°. SER occurs only at the upper and lower tips of 
the constant energy band contour when \ky\ is large. It should be 
emphasized that valley-selective reflection is only possible when elec- 
trons undergo retro -reflection. The specularly reflected electrons are 
always divergently reflected and the valley states are mixed since the 
upper and lower 'tips' of the constant energy band contour at large \ky\ 
are approximately identical for both K and K' valleys. In Fig. 6d, the 
RER trajectory of electrons injected from a point source is shown. The 



SCIENTIFIC report; | 2 : 1013 | DOI: 1 0.1 038/srepOl 01 3 



4 




-2-1012 _2 -1 0 1 

kx/h) kx/ko 




Figure 6 | The electron excitation spectrum of BLG/S at (a) a IC valley and (b) a K' valley with £p = 0.5£o and Ej^ = £q. The RER constant energy band 
contour of the iC valley (red dashed curve) is significantly 'smoother' than that of the K' valley (green dashed curve). At large I ky\ (SER regime), Kand K' 
band contours become approximately identical; (c) Reflection angles of iC(red curve) and iC' (green curve) electrons. electrons are predominantly 
focused via smaller angle (< 20°) than that of the K' electrons (—80°). (d) RER trajectory of electrons emitted from a point source situated at P (denoted 
by yellow triangle). Blue rays represent the incident electrons. The interface acts as a valley-selective dual-focus electron mirror with iC electrons (red rays) 
being focused further away from the interface than the K' electrons (green rays). 



BLG/S junction acts as a valley- selective dual focusing electron mirror 
with K electrons being focused further away from the junction due to 
the 'smoother' constant energy band contour of K reflection states 
while K' electrons are focused closer to the junction due to the 
'sharper' constant energy band contour of K' reflection states. 

Discussion 

In Fig. 2, we show that the anisotropic constant energy band contour 
of BLG allows RER to occur. The reversal of Vy group velocities upon 
reflection, as shown in Fig. 3, further illuminates the peculiar RER 
phenomena in BLG. It should be emphasized that the RER is more 
profound in the low energy regime of ~ Sq- This requires high 
quality samples with carrier concentrations n ~ 10^^ cm"^ ^ which is 
experimentally achievable^^. RER does not survive in energy regimes 
higher than > 10.9eo since the reversal of Vy upon reflection can no 
longer occur as the anisotropy of the constant energy band contour is 
lost. 

The absence of AR in a BLG/S junction at low energy regimes is 
in agreement with previous study on BLG/S junction utilizing a 
four-band model without low energy band anisotropy^^ It was 
shown that AR in BLG is a small effect in the order of Ek/yi. 
Here we further demonstrate that AR is not only a small effect 
but is completely absent in the low energy regime. The absence of 
AR is related to the Klein tunneling effect and Berry phase in 
BLG. In BLG with isotropic and parabolic energy dispersion, elec- 
trons incident normal to a potential barrier are backscattered with 
unit-probability since the transmission mode normal to the inter- 
face is evanescent as encoded in the 2n Berry phase of its spinor 
wavefunction^'*'^^. The 2n Berry phase is also retained when the 
low energy anisotropy is included^^. Due to the large Fermi mis- 
match at the junction (£|f'*^e|f^^'*), the transmission modes are 
always perpendicular to the interface for any incident angles and 
hence instead of perfect electron backscattering at only normal 
incidence, electrons reflect with unit probability at any incident 
angles in a BLG/S junction. Mathematically, this can be explained 



by the unique spinor-form of the BLG low energy electron wave- 
functions. The transmitted spinor-wavefunction in the supercon- 
ductor side for electron-like and hole-like excitations are, 
respectively, 

(3) 

(i r) (i r) 

where u and v are the BCS coherence factors, <Aa| » ybI corre- 
sponds to the wavefunction amplitude at sub -lattice A | and B | 
respectively. The superscript (/, r) denotes incidence and reflection 
states respectively. Upon transmission, the normal-transmitted 

spinor follows the equality = This essentially uncouples 
the electron- reflection coefficient from the other transmission/ 
hole- reflection coefficients, resulting in our angular-independent 
result: Re = I (see Supplementary Information). This aspect is in 
stark contrast to the single layer graphene/superconductor junc- 
tion (SLG/S) where electrons reflect in an entirely opposite way: 
AR occurs with unit-probability only at normal incidence^. The 
perfect normal-transmission of electrons in graphene is a sig- 
nature of the 71 Berry phase of the massless Dirac fermions^^'^^. 
The superconducting wavefunction in SLG/S is in the form of 

(/ r) (i r) 

rA ' rB where subscripts A and B correspond to the wave- 
function amplitudes at the two inequivalent sub -lattices A and 
B respectively. Upon transmission, the spinor follows a rather 

different relation of xj/^^ = — ij/^g^ . The hole and electron reflection 
coefficients are mathematically entangled with each others in this 
case and hence the AR reflection coefficient is angularly depend- 
ent on the electron incident angle. 

In conclusion, graphene has again surprised us by allowing the 
peculiar RER to occur in its bi-layered form. The occurrence of RER 
in BLG while avoiding the violation of the momentum conservation 
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rule represents a conceptually simple yet astonishing example high- 
lighting one major difference between classical mechanics and 
quantum mechanics. In a BLG/S hybrid junction, the In Berry phase 
of low energy electrons in BLG manifests itself in a striking manner: 
AR is completely absent regardless the angle of incidence. This makes 
BLG/S a rather ideal platform to demonstrate RER since the com- 
peting AR process is fully eliminated. The absence of AR in BLG/S 
junction can be experimentally verified as the absence of supercur- 
rent flow in the low Fermi level regime. Finally, the valley- inequiva- 
lent RER phenomena in BLG/S is analogous to a valley- selective 
Veselago electron-mirror. The valley- selective dual-focusing beha- 
vior suggests a potential application in the field of valleytronics. 

Method 

See Supplementary Material for derivation of the reflection coefficients. 
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